set term png
set output "mc_markov_demon_E.png"
muB = 2
f(kT) = (-muB*exp(muB/kT) + muB*exp(-muB/kT))/(exp(muB/kT) + exp(-muB/kT))
set xlabel "kT"
set ylabel "<E>"
plot [0:50] f(x)

set output "mc_markov_demon_Sz.png"
g(kT) = (exp(muB/kT) - exp(-muB/kT))/(exp(muB/kT) + exp(-muB/kT))
set xlabel "kT"
set ylabel "<Sz>"
plot [0:50] g(x)

